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We have done a study with small, imperfect Hubbard rings with exact diagonalization. The results 
for few-electron rings show, that the imperfection, whether localized or not, nearly always decrease, 
but can also increase the persistent current, depending on the character of the imperfection and the 
on-site interaction. The calculations are generally in agreement with more specialized studies. In 
most cases the electron spin plays an important role. 



PACS numbers: PACS: 73.23.Ra,73.21.Hb,73.21.La 



I. INTRODUCTION 

The regime of experimental studies in small semicon- 
ductor heterostructures has gone from mesoscopic to 
nanoscopic within a few years. Especially electron traps 
have been under intensive study because of the discrete 
nature of energy levels and the resulting analogy with 
real atoms. The trapping potential has varied between 
parabolic (quantum dot) 1] and ring-like (quantum ring, 
QR) 1^ 1^ 0| ' While quantum dots fascinate because of 
the similarity with real atoms, rings combine this simi- 
larity with an always captivating ring-like geometry. 

Rings are often studied with a relation to the persis- 
tent current (PC) 0, an equilibrium current that arises 
when an Aharonov-Bohm flux is piercing the ring 0. 
Few-electron nanoscopic rings introduce many-body ef- 
fects not observable in mesoscopic rings, such as the 
fractional 0o /iV-periodicity of the persistent current . 
These systems have been studied theoretically, both in 
the single-particle B, Q as well as in the many-body pic- 
ture 0, lo. [T3L Il4| . With continuum models, partic- 
ularly in the analytical approach, it is laborious to in- 
troduce imperfections to a perfect ring, especially when 
non-perturbative treatment is required. 

In addition to continuum models, lattice models have 
been applied to QRs 0, 0, E, 0| ■ If interactions are 
taken into account, the Hubbard model |23| is probably 
the most investigated model. Now, apart from being a 
toy model of mathematical physics |2lj . the purpose of 
this paper is to show, that a slightly generalized Hubbard 
model can also be an extremely valuable tool in getting 
the first crude idea of what is going on in small, numeri- 
cally exactly diagonalizable systems. With the Hubbard 
model and exact diagonalization it is easy to introduce 
the effects, including for example an impurity, disorder, 
magnetic- and electric fields and external leads, that are 
much harder to include with other models. Here one 
does not need to limit to a certain region in the param- 
eter space but can explore all the possible values of the 
parameters and we get an illustrative representation of 
the results. 

The calculations show, that imperfections mainly de- 
creases the persistent current. The interaction, however, 
can introduce coupling of the two spin-currents and re- 



sult in increase of the PC in certain situations, especially 
when the imperfection is of localizing nature. We con- 
clude that spinless models cannot show these effects. 



II. THE MODEL 

We use the Hubbard Hamiltonian with a pure 
Aharonov-Bohm flux without a Zeeman term. In the 
presence of a vector potential the hopping integral is 
modified by a phase factor and the Hamiltonian reads 



+ U J2t fltinil +J2ia 



(1) 



where we have generalized for site-dependent hopping in- 
tegral as well as for one-body on-site energies Ci . Aij rep- 
resents the portion of the Aharonov-Bohm flux covered 
by the hopping ? j, so that we could write the phase 
factor as exp[—i(2n / L){(f)/ (po)], where L is the number of 
sites in the ring, (jj is the magnetic flux and (j)o — h/e the 
flux quantum. Energy scale is fixed by setting all tij — 1 
unless otherwise stated. Figure la shows examples of 
Hubbard rings considered in this research. 

The dimension of H is (^^) (^^) , where N^y is the num- 
ber of spin-tT -electrons, growing strongly with L and N^, 
limiting the system size. The diagonalization is done with 
standard library routines and the numerics do not involve 
severe difficulties. 

The persistent current of an eigenstate ipm in a sin- 
gle ring is usually viewed via the expression Im{(f>) — 
-dE„Jd(j) m ig. In the presence of multiple rings 
and different currents, however, we have to use the cur- 
rent operator. From the relation above and from the 
Feynnian-Hellmann theorem we obtain [27| 



Jklo 



Im[e-^(2VL)(0/0o)ct^c^^] 



(2) 



for the spin-fj -electron current operator between the (ir- 
relevant) sites k and /. The persistent current with the 
definition = —dE/d<j) is the sum of the two differ- 
ent spin currents. Using this operator is beneficial also 
because we obtain the current directly by taking the ex- 
pectation value of the operator jSJ in a given state. 
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FIG. 1: (a) Structures of some of the Hubbard rings used 
in the calculations. The size of the point represents the site 
energy. The dashed line represents a weaker (or stronger) 
hopping propability. (b) Persistent current of a perfect ring 
with L = 7 and A'^j = A*'; = 2 as a function of the interaction 
U and the flux (p. Note the periodicity change from <j!>o to 
4>o/N^ — </>o/4 as the interaction increases. 



The persistent current is a periodic function of 4> with 
periodicities (j)o, </>o/2 and (po/Nf,, increasing with the in- 
teraction, as shown in Fig . lb. The periodicity has been 
studied before 0,0,|23|> the purpose of Fig. lb is 
just to show that it makes sense to characterize the per- 
sistent current by fixing (f> as long as it does not contain 
any discontinuities. Like many other authors Il3|. we 
have used the value cj) = O.250o and the ground state 
persistent current throughout the paper. 

It is suggestive, that in the Hubbard model the quasi- 
one-dimensionality can be mimicked by adjusting the in- 
teraction term: A Hubbard ring with multiple channels 
(Fig. la) and J7 ^ 1 is mimicked by a single channel and 
a relatively small U ■ Thus the Hubbard U in these 
ring-like structures can be read either as an interaction 
as such or as a measure of the one-dimensionality of the 



ring. 



III. RESULTS AND DISCUSSION 

The effect of a single impurity on PC has been studied 
by several authors [ilia, H^. Fig.2a shows PC influ- 
enced by an on-site energy ei and the interaction. With 
zero impurity, the current decreases when the interac- 
tion U (or one-dimensionality) increases. This happens 
because with increasing U the Hamiltonian prevents 
electrons moving independently and eventually only ro- 
tation as a rigid rotor is possible. These observations are 
in agreement with more specialized papers |T3i il6.j . 

The current always decreases monotonously with in- 
creasing positive ei. Investigating at the Fig. 2a more 
carefully with a fixed positive impurity, one can see that 
there is a slight increase of the current with the interac- 
tion. This result, confirming earlier calculations |3^l35l|. 
tells that a greater interaction makes the electron system 
more correlated, pushing the impurity-affected electrons 
more effectively. Because there are observations contra- 
dicting with this result saying that in the presence 
of impurity, the interactions suppress PC even further, 
we hasten to add that with ei > the current with large 
enough interaction is slightly smaller than the current 
without interaction. Note that in the case of, for ex- 
ample, the Coulomb potential, the interaction is either 
on or off, and this kind of gradual increase in the inter- 
action strength is usually not considered. However, as 
discussed above, the increase of electron-electron inter- 
action corresponds to decrease of the width of the quasi- 
one-dimensional ring. 

With often disregarded attractive impurity and fixed 
U one can see a local maximum of / with specific ei. 
Here the impurity localizes electrons, and if the depth of 
the impurity potential is equal to the interaction, U + 
ei = 0, the electrons with the opposite spin as an effect 
do not feel any extra on-site energy and thus pass the 
impurity unaffected. Either increasing or decreasing U 
(or ei) makes the effective on-site energy repulsive or 
attractive, decreasing /. Away from this region of local 
maximum (especially with [7 = 0) the attractive impurity 
decreases the current more effectively than the repulsive 
one. 

A similar effect is seen in Fig. 2b, where the "impurity" 
is now a different hopping integral ti2. The current goes 
to zero with ti2 as it should, but surprisingly, it has a 
maximum near ti2 ~ 1, above which it decreases. The 
physics is the same as with attractive impurity, because 
the large negative kinetic energy of the strong link local- 
izes an electron pair to the corresponding sites, creating 
a similar kind of blocking effect; notice that the form of 
Fig. 2b with ti2 > 1 is the same as in Fig. 2a with ei < 0. 

Because we have an equal number of electrons with the 
opposite spin, the different spin-currents, by symmetry, 
are the same. But if N-^ ^ N^^, they are in the opposite 
direction. This can be seen in Figs. 2c and 2d, which show 
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FIG. 2: Effect of a single impurity on the persistent current at a fixed value of flux {(f> = 0.2500)- In (a) the impurity is 
described with a onsite energy ei, in (b) with a hopping parameter ti2 {L = 7, Nt^ = = 2). (c) and (d) show separately the 
spin-up and spin-down currents in the case of three electrons (L = 7, A^'i = 2, = 1) and a single impurity described with ei. 



spin-up and spin-down-currents as a function of U and 
the single impurity strength for = 2 and Ni — 1. The 
opposite signs and relative magnitudes of PC with zero 
interaction and impurity, i.e. /|(0, 0)//j^(0, 0), can in fact 
be easily explained also with non-interacting electrons in 
a continuous, strictly one-dimensional ring. The effect of 
interaction is to 'grab' the spin-down electron to move 
to the same direction as the spin-up electrons, leading 
eventually to the rigid rotation of the whole electron sys- 
tem. However, to the total current (/| +Ii) the impurity 
has the same effect in the three-electron case as in the 
four-electron case. With other odd numbers of electrons 
the situation is similar to that of the three-electron case. 

A random external potential at the ring is often called 
disorder and its effect on PC and the energy spectra is 
largely studied subject mainly due to the relevance in ex- 
perimental samples pA ll7lll8l |. Fig. 3a shows an ensemble 
averaged I{U,W), where W is the measure of the dis- 
order, defined as W = max{{€i}) — min{{€i}), and {ei} 
are random on-site energies (the enesemble consisted of 
60 rings). The current decreases monotonously with W, 
but the role of the interactions is again not trivial, since 
a given disorder strength gives the current a maximum 



with U ^W. The increase of PC with [/(< 1) has been 
confirmed by other authors [33. l35l| . Physically the max- 
imum can again be explained by the effective potential- 
smoothening that originates from the disorder-localized 
electrons and their repulsive interaction towards the elec- 
trons of opposite spin. The impact of interactions can be 
different in spinless models [13 (where inter-site interac- 
tion is required, e.g. a Coulomb-like interaction term 
0). But at this point we want to stress that if the spin- 
currents, potentially even of opposite sign, are present, 
then the effect of interaction, by coupling the two cur- 
rents as described above, is quite subtle and differs from 
the effect in spinless models. 

By comparing carefully Figs. 3a and 3b, the resem- 
blance is obvious. Fig. 3b shows the same plot as Fig. 3a 
but now with disorder strength W replaced by the 
strength Ve of an electric field parallel to the plane of 
the quantum ring. The physics remain the same: while 
the field 'inclines' the ring so that some electrons tend 
to roll down to the other end of the ring, the interac- 
tion smoothens the way to the electrons of the opposite 
spin, creating an optimum interaction for given Ve- On 
the other hand, for a given U an increasing electric field 
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FIG. 3: Effect of random on-site energies (a) and electric 
filed (b) on the persistent current ai tf) — O.250o {L = 7, 
iV| = A'^l = 2). W is a measure of the disorder of the onsite 
energies and Ve is the strength of the electric field parallel to 
the plane of the ring. 

destroys the Aharonov-Bohm oscillations, since the wave- 
function does not extent throughout the whole ring |37j . 

Finite temperature behaviour of the persistent current 
is more intuitive, as seen in Fig. 4a. The current de- 
creases monotonously with U at all temperatures, and 
with a fixed U the current is almost constant for small 
enough T, until it decreases essentially to zero in the 
temperature range that is given by the energy gap be- 
tween the ground state and the excited states. This is 
because the directions of the currents of the exited states 
are frequently opposite to the current of the ground state. 
The observations agree with the non-interacting picture 
of Refs.[lll3. 

As is the temperature, so are the external leads an 
inevitable imperfection in externally tuneable quantum 
rings [^Q. With a lattice model one could imagine mod- 
elling these with an attractive on-site energy, larger hop- 
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FIG. 4: (a) Temperature dependence of the persistent current 
in a perfect ring with with L = 7 and = A''| = 2 {ksT 
measured in units of t). (b) The energy spectrum Ei{(j)) of 
a ring with six sites and a stub with two sites (in Fig. la is 
shown a related system with a ring of seven sites and a stub 
of three sites). The system has four electrons {N-[ = Ni —2) 
and U = 1000. The low-energy spectrum is nearly identical 
with a spectrum of three electrons in a pure ring without the 
stub. 



ping integral or maybe with a smaller interaction term 
that is due to the larger space available in the vicin- 
ity of the leads. The ring-stub -system, a schematic of 
which is shown in Fig, la, has been studied by contin- 
uum models I2E 1^- It was found, that the stub 
can create standing waves which are not affected by the 
magnetic flux piercing the ring. With our approach, in- 
cluding now also the many-body effects, it is found that 
with large enough interaction, the stub of definite length 
can localize an integer number Ngtub of electrons, and as 
a result the remaining electrons in the ring exhibit nor- 
mal AB-oscillations with Ne — Nstub electrons, as if the 
stub would be absent. This is shown in Fig. 4b, where we 
have a system of four electrons, but the low-energy spec- 
trum is essentially the same as the pure-ring spectrum 
with just three electrons. If the stub-length would be in- 
creased by one site, the standing wave would not have a 
node anymore at the stub-ring interface and the localized 
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electron would start interfering with the electrons in the 
ring, suppressing the current. 

If the stub is connected to another quantum ring, as 
depicted in Fig. la, we get coupled quantum rings, which 
could be experimentally realistic in a dense ensemble 
of self-assembled quantum rings 2). Studied with non- 
interacting continuum electrons _33], the relative direc- 
tions of PC in coupled rings of different radius were found 
to depend on the magnetic flux. If we define a reference 
ring with length to enclose a flux (j}^, a ring with 
length L then encloses a flux (p = {L/Lr)'^(l)r, and we can 
have rings of different sizes by applying the correspond- 
ing fluxes to the phase factor of the hopping integral in 
the Hamiltonian It is found that the currents in the 
rings of length Li and L2 indeed run in the opposite di- 
rection if e.g. the flux through the first ring is less and 
the flux through the second ring is more than integer 
multiple of a flux quantum |4l|. Though we have two 
independent currents, it turns out that the derivative of 
the total energy resembles the sum of the currents very 
closely, and thus may lead to large periodicities if L1/L2 
differs only slightly from unity. 

Most of the figures show results for rings with L — 7 
and Ne = 4, but it is important to point out that also 
other rings were investigated, giving similar results. We 
want to defend this choice by the notion that in the limit 
of large U, the Hubbard model becomes the Heisenberg 
model with an effective coupling constant Je/ / [iS ' 
which, for large L, scales as Je// oc L~^. This, together 
with our calculations, shows that the energy separations 
-and consequently the persistent current- decreases with 
L. This implies that the general trends depicted in the 
figures tend to diminish as the number of empty lattice 
sites is increased. Furthermore, these few-electrons sys- 
tems, because of the local interaction, the model does not 
necessarily compare to real continuum models better if L 
is increased. 



IV. CONCLUSIONS 

To conclude, we have studied small, imperfect Hub- 
bard rings with exact diagonalization. The results show, 
that the imperfection, whether localized or not, almost 
always decrease the persistent current. The interaction, 
however, can introduce coupling of the two spin-currents, 
make the effective potential felt by opposite-spin elec- 
trons smoother and result in increase of PC in certain 
situations, especially when the imperfection is of local- 
izing nature. We note that these effects cannot be seen 
with spinless electrons. The decrease in PC as a function 
of impurity potential starts with zero slope. The cur- 
rent decreases with temperature monotonously in a scale 
given by the energy gap above the ground state. 

By comparing with other studies, it is shown that even 
these very small rings give the same physics and phenom- 
ena that is obtained with more involved Hubbard-model 
calculations, including e.g. Bethe-Ansatz -calculations. 
Furthermore, many realistic investigations with contin- 
uum models leave the phenomena obtained with this ap- 
proach mainly intact. 
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